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We present searches for charged Higgs production in decays of top quarks and also pair production of doubly 
charged Higgs boson decaying to di-tau, di-muon, and muon + tau final states. The searches are performed in 
proton-antiproton collisions at */s = 1.96 TeV using an integrated luminosity of up to 7 fb -1 collected by the 
CDF and DO experiments at the Fermilab Tevatron Collider. We find no evidence for charged Higgs production 
and set limits on the production cross-section for a variety of theoretical models. This represents the first search 
for pair production of doubly-charged Higgs bosons decaying into tau leptons at a hadron collider. 



1. Introduction 

Charged and doubly charged Higgs appear in various extensions to the Standard Model, SM. In order to 
generate a charged Higgs one needs to have two Higgs doublets, describing the Higgs sector, compared to 
the one in the SM. One set of theories of particular interest are supersymmetric, SUSY, theories. These are 
extensions to the SM which predicts a new symmetry between bosons and fermions. This theory has several 
advantages over the SM such as the introduction of a dark matter candidate, a solution to the hierarchy problem 
and the potential for GUT scale unification. In its simplest form it is the MSSM, the minimal supersymmetric 
standard model, which after electroweak symmetry breaking predicts 5 Higgs boson, 3 neutral and 2 charged. 
Charged Higgs also appear in extended SUSY models such as the NMSSM, the next to minimal supersymmetric 
standard model [2J. At tree level the Higgs sector of the MSSM can be described by two parameters which are 
chosen to be, tan/3, the ratio of the vacuum expectation of the two Higgs doublets and, Ma, the mass of the 
pseudo-scalar Higgs boson. Beyond tree level radiative corrections bring in dependance on more than just Ma 
and tan/3 so limits are set at specific benchmark regimes in the Myt-tan/3 plane. To produce a double charged 
Higgs the SM must be extended further with the inclusion of a Higgs triplet, this occurs in models such as the 
Left Right Symmetry Model [|[ and the Little Higgs models 

2. Charged Higgs 

The decays of charged Higgs depends on both the model considered and the mass range predicted. Various 
scenarios have been studied at both DO and CDF, for non-extended two Higgs doublet models, such as the 
MSSM, the charged Higgs is expected to couple strongest to the top quark. Therefore analysis studied can be 
separated in cases where the charged Higgs is lighter than the top quark and where it is more massive. 

2.1. Two Higgs doublet models: M top > M H 

For the case where the mass of the charged Higgs is smaller than the mass of the top quark, the production 
of a charged Higgs at the Tevatron will be via t —> H + b. The charged Higgs can either decay leptonically, 
hadronically or to a mixed hadron lepton final state. This additional decay of the top quark will compete with 
the SM process t — » W + b. Due to the different decays of the W + and H + a search for the charged Higgs can 
be preformed by looking at a difference in the distributions of events between the final states as compared to 
what is predicted by the standard model. Searches had been performed for a leptophobic model where it is 
assumed that the BR(iJ + — > cs) ~ 1. This was preformed at CDF with 2.2 fb _1 of integrated luminosity and 
at DO with 1.0 fb _1 DO has also conducted a search for both the leptophilic case where it is assumed that 

BR(H + tv) ss 1 and a strangephilic model where BR(H + cs) + BR(H + ->• rv) m 1 @ with 1.0 fb _1 of 
integrated luminosity. 

2.2. Two Higgs doublet models: M H > M top 

If it is assumed that the mass of the charged Higgs is greater than the mass of the top quark, then the charged 
Higgs will be produced through pp annihilation and will decay to a tb pair. This is the same final state as the 
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Figure 1: 95% confidence limits on the cross section x branching ratio of a charged MSSM Higgs boson shown in the 
M H + -tan/3 plane for a strangephilic model of charged Higgs coupling. The expected limit is shown by the red line and 
the hatched blue region is the region excluded at the 95% confidence level. 

standard model process W + — > tb therefore the rate of single top production at the Tevatron compared to 
that predicted by the SM would be effected and a bump would be expected in the invariant mass distribution 
of the tb final state. DO performed a search for a charged Higgs boson is this region between the masses of 
180 < M H + < 300 GeV using 0.9 fb -1 of integrated luminosity. No evidence of a charged Higgs boson was 
found and limits where set on a(pp —> H + )x BR(H + —> tb) for three different types of two Higgs doublet 
model. In Type One two Higgs doublet models only one of these doublets couples to fermions and a region in 
the plane of the charged Higgs mass and tan/3 and be excluded [7|. This can be seen in Figure 2. 
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Figure 2: 95% confidence limits on the cross section x branching ratio of a charged MSSM Higgs boson shown in the 
M H + -tan/3 plane for a heavy charged Higgs decaying into tb pairs. The region excluded is shown in green. 



2.3. Extended two Higgs doublet models 

Charged Higgs boson also appear in extended Higgs models such as the NMSSM, this model extends the 
MSSM by the addition of a singlet Higgs boson on top on the two doublets in the MSSM. After electroweak 
symmetry breaking this results in 7 Higgs bosons, 5 neutral and 2 charged. The NMSSM has several advantages 
over the MSSM, it provides a natural scale for the /x parameter, supersymmetric mass parameter, it alleviates 
the little hierarchy problem and for regions of the NMSSM parameter space where the mass the lightest pseudo- 
scalar Higgs boson, a, is less than twice the mass of the b quark, m a < 2m;, then the LEP limits can be 
avoided Q. Within the NMSSM, for large regions of the parameter space Higgs-to- Higgs decays dominate over 
decays to other particles. A search has been performed at CDF for a NMSSM charged Higgs boson with 2.7 
fb -1 of integrated luminosity Q. It was assumed that the charged Higgs will be produced from the decay of 
a top quark and will decay into a light psuedo-scalar Higgs boson. The mass of scalar Higgs is assumed to be 
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less than twice the mass of the b quark and hence it will to decay to a rr pair, t — > H + b — » W ± ab — > W^rrb. 
Limits are determined on the t — » H + b branching ratio for various H + and a masses by preforming a fit to the 
isolated track pr spectrum. The determined limits on the BR(t — > H^b) as a function of the charged Higgs 
mass are shown in Figure 3. 
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Figure 3: 95% confidence limits on the BR(t —> H b) as a function of the charged Higgs mass. Four different values of 
the pseudo-scalar Higgs mass are shown. 



3. Doubly charged Higgs 

Doubly charged Higgs appear in extensions to the SM in which introduce a Higgs triplet. After electroweak 
symmetry breaking this results in Higgs boson with a double charge. Higgs triplet models that have been 
covered by analyses at the Tevatron included the Left Right Symmetry model, which predicts a new symmetry 
between left and right handed particles Q • In this model there is both a left-handed and a right-handed doubly 
charged Higgs boson, H^,!!^ where is predicted that the production cross section of the right handed Higgs 
is approximately half that to the left-handed state due to the difference in the coupling to the Z boson where 
the left- and right-handedness of the doubly charged Higgs boson refers to the handedness of the particles that 
couples to and not of the Higgs itself. Higgs triplets also appear in Little Higgs models, with a type 2 see-saw 
mechanism where they are used as a production mechanism for neutrino masses 4] . For a normal hierarchy of 
neutrino masses below approximately 10 meV, these models predict that the doubly charged Higgs will have 
approximately equal branching ratios to /x/Lt, rr and fir final states [9j. In (3 — 3 — 1) gauge symmetric models 
which predict additional heavy exotic quarks and leptons that provide anomaly cancellations, double charged 
Higgs that decays dominantly to taus [l(|. DO has performed a search looking for the pair preproduction of 
iJ ±=t decaying into r ± r ± ,/x ± /z ± and /i^r 1 final states using up to 7.0 fb _1 of integrated luminosity. Detailed 
description of this analysis can be found in [ll[ . This is the first search for doubly charged Higgs decays to a pure 
hadronic tau final state, i/ ±± — > r ± r ± , at a hadron collider. Limits were set for both model dependent scenarios, 
i.e. assuming a equal BR to /i/i, tt and \it final states as predicted in @ and for model independent ones. These 
were BR(# ±:b -> t ± t ± ) = 1, BR(i? ±± -> ^ ± A* ± ) = 1, BR(iJ ±i: I^t*) = 1 and BR(i? ±± r ± T ± ) + 
BR(i/ ±± — > /i ± /i ± ) = 1, where BR(i? ±± — > r ± r ± ) is ranged over 10 different value from 100% BR to taus 
to 100% BR to muons. All the model independent limits are determined for both and right- and left-handed 
Higgs. 

3.1. Analysis Summary 

This analysis starts with the identification of muons and hadronic taus are identified. For a reconstructed 
muon, hits in the muon chambers are matched to tracks in the central tracking chambers, these are then required 
to be isolated in both the calorimeters and the central tracking detectors. Hadronic tau decays are first split 
into three types based on the energy clusters in the calorimeters and the number of tracks. Type 1 and 2 are 
1-prong decays with energy deposited in the hadronic calorimeter (type 1) or in both the electromagnetic and 
hadronic calorimeters (type 2). Type 3 taus are 3-prong decays with an invariant mass below 1.7 GeV and 
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energy deposits in the calorimeters. For each of the different types a Neural Network, NN, is trained using 
Z — y tt decays as signal and multijet events predicted from data as background, to separate hadronic tau 
decays from jets. 

Events are required to have two hadronic taus, 77,, and at least one isolated muon. The following selection 
criteria must be satisfied; the pseudo-rapidity of selected muons and taus must be |?/'| <1.6, \n Tl < 2 \ < 1.5, 
additional taus must have |77 r | < 2. The transverse momentum of the objects must satify pt^ > 15 GeV and 
pTj, > 12-5. All Th and muons are required to be separated by AR > 0.5, where AR — (A</>) 2 + (An) 2 , for 
the two leading px taus then AR > 0.7. The sum of the charges of the highest pr muon and the two highest 
Pt taus is required to be Q = 53i=u Tl t 2 * = as * s ex P ec ted for signal. After application of the described 
cuts, the main backgrounds contributions are from diboson and Z — ¥ tt processes, with smaller contributions 
arising from Z —} fifi, W+jets, tt and Z — > ee. All background process are simulated with ALPGEN [l2j with 
showering and hadronization provided by PYTHIA fl3l | except for diboson and signal events which a re g enerated 



by PYTHIA. Tau decays are simulated using TAUOLA which correctly model tau polarization 14J. For all 
MC samples GEANT [15| is used to correct for detector effects. The contribution due to multijet background 
is determined to be negligible. 

To enhance the discrimination of the signal sample over the background the selected events are split into 
four different samples. These are based on the number of tau and muons in the final state and the charged 
correlation between the the two taus. For two samples it is required that there is exactly two taus and one 
muon and these are split further depending on whether the charges of the two taus are the same or opposite. 
The other two samples require exactly 3 taus and 1 muon or exactly 2 taus and 2 muons. The invariant mass 
of the two highest pt candidates for the two samples with exactly 2 taus and 1 muon are shown in Figure 4. 
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Figure 4: The invariant mass of the two highest pr taus for the selection with 2 taus and 1 muon for (a) q(ri) = q(r2) (b) 
Q(Ti) = — q{T2)- The crosses show the data, the stacked histogram the simulated background and the solid lines expected 
signal for three values of the H ++ branching ratio and a mass of 120 GeV. 



The four different samples are dominated by different background processes hence splitting improves the 
sensitivity to the signal. The like charged sample is dominated by Z/j* events, mainly —¥ tt with 

contributions from Z/j* — > /i/i and VF+jets. The opposite charged sample is dominated by diboson events 
mainly WZ — > /ui/e+e - with electrons misidentified as taus. Smaller contributions arise from Z/j* and W^+jets 
events. Table I shows the expected number of signal and background events for the four different signal samples. 

As no significant excess in the data over the background is seen, limits were set on the H ++ production cross 
section using a modified frequentist approach [16]. The invariant mass of the two highest pt taus, M(ti,T2), is 
used to discriminate signal from background. 

The dominant sources of systematic uncertainties arise from, the uncertainty on the measured integrated 
luminosity, the uncertainties on the muon and tau identification including the uncertainty for applying the NN r 
and the uncertainties on the background cross sections for Z/j*, W^+jets, tt and diboson production. There 
are also uncertainties from the trigger efficiency and on the signal acceptance due to the parton distribution 
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Table I: Numbers of events in data, predicted background, and expected signal for M(H^ ) = 120 GeV, assuming the 
NLO calculation of the signal cross section for BR(//^ ± -> r ± r ± ) = 1, BR(i?* ± -> M ±r± ) = 1, and BR(i/* ± -> 
r ± r ± ) = BR(i/^ :± — > fj, fj, ) — BR(/f^ :± — > /n r ) = 1/3. The numbers are shown for the four samples separately, 
together with their total uncertainties. 
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Figure 5: The upper limit on the H^^H^^ pair production cross section. For the case where BR(if ±± rr) = 
BR(// ±± — ► fifi) = BR(i/ ±± — > t/jl) — 1/3. The black solid and dashed lines show the observed and expected limits 
respectively. The yellow and green regions are the one and two standard deviation bands and the blue band is the 
predicted cross section with associated uncertainty. 



function. 

Limits are set assuming several different values of the H ++ branching ratio (a) BR(i7 ±± — > r ± r ± ) = 1, 
(b) BR(i7 ±± -> r^) = 1, (c) BR(i/ ±± t^t*) = BR(i? ±± r ± /* ± ) = BR(i7 ±± -> ^ ± A* ± ) = 1/3 
and (d) assuming BR(i7 ±± -> r ± r ± ) ± BR(F ±± ^ ± m ± ) = 1. For (d) the BR(i? ±± -> r ± r ± ) is ranged 
in 10% increments from 100% BR to taus to 100% BR to muons. In order to cover all this parameter space 
this analysis is combined with the a search from the DO collaboration pjj for the pair production of H ++ 
decaying into 4/x final state, using 1.1 fb _1 of integrated luminosity. The distribution of the invariant mass of 
the two highest pr muons, as determined by this analysis, including all systematics and their uncertainties and 
correlations here included in the limit setting procedure. For all the above scenarios, limits were set for both 
a left-handed and right-handed doubly charged Higgs apart from scenario (c) as the motivation is a model in 
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which only left-handed Higgs exist. Figure 5 shows the determined limits for scenario (c). Figure 6 shows the 
limits for scenario (d) for both a left-handed and right-handed doubly charged Higgs. The limits determined 
are summarized in Table II and are the most stringent limits on a pair produced doubly charged Higgs in these 
decay channels and the only limits set a hadron collider in the case when it is assumed BK(H ±:t — s- r ± r ± ) = 1. 
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Figure 6: The expected and observed exclusion regions at 95% confidence level in the BR(ff J 
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Higgs. The yellow shows the uncertainty on the NLO calculation of the signal cross section. 



Table II: Expected and observed limits on M (// ±± ) (in GeV) for left and right-handed H ±i: bosons. Only left-handed 
states are considered for the model that assumes equality of branching fractions into tt, (j,t, and fj,fj, final states. We 
only derive limits if the expected limit on M{H ) is > 90 GeV. 
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